The Optical/Near-Infrared Light Curves of SN 2002ap for the 

First 1.5 Years after Discovery 
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ABSTRACT 

Late-time BVRIJHK photometry of the peculiar Type Ic SN 2002ap, taken between 2002 
June 12 and 2003 August 29 with the MAGNUM telescope, is presented. The light curve decline 
rate is derived in each band and the color evolution is studied through comparison with nebular 
spectra and with SN 1998bw. Using the photometry, the OIR bolometric light curve is built, 
extending from before light maximum to day 580 after explosion. The light curve has a late-time 
shape strikingly similar to that of the hypernova SN 1998bw. The decline rate changes from 0.018 
mag/day between day 130 and 230 to 0.014 mag/day between day 270 and 580. To reproduce 
the late-time light curve, a dense core must be added to the 1-D hypernova model that best fits 
the early-time observations, bringing the ejecta mass from 2.5 Mq to 3 Mq without much change 
in the kinetic energy, which is 4 x 10 51 ergs. This is similar to the case of other hypernovae and 
suggests asymmetry. A large ff-band bump developed in the spectral energy distribution after 
~ day 300, probably caused by strong [Si I] 1.646/xm and 1.608^m emissions. The near-infrared 
flux contribution increased simultaneously from < 30% to > 50% at day 580. The near-infrared 
light curves were compared with those of other Type Ib/c supernovae, among which SN 19831 
seems similar to SN 2002ap both in the near-infrared and in the optical. 

Subject headings: supernovae: general — supernovae: individual (SN 2002ap) — supernovae: photometry 
1. Introduction 
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There is growing interest in Type Ib/c super- 
novae (SNe) following the discovery of a subclass 
that display very broad-lined spectra, indicating 
the existence of ejecta expanding with velocity 
> 0.1c. Models of the optical spectra and the 
light curves (LCs) of this high-velocity subclass, 
also called hypernovae (HNe, e.g., Nomoto et al. 
2004), have concluded that these are the explo- 
sions of massive C+O stars, producing up to 10 
times the kinetic energy of normal core-collapse 
SNe (e.g., Iwamoto et al. 1998, 2000; Mazzali et 
al. 2003; Deng et al. 2005). It is unclear why 
some HNe are apparently associated with GRBs, 
i.e., SN 1998bw and GRB 980425, (Galama et al. 
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1998), SN 2003dh and GRB 030329 (Stanek et al. 
2003), and SN 20031w and GRB 031203 (Malesani 
et al. 2004), while others are not. 

The Type Ic SN 2002ap, discovered in M74 
(Nakano et al. 2002), is one of the nearest SNe 
in recent years. Its host galaxy has a distance 
modulus of only 29.5+0.2 (Sharina, Karachentsev, 
& Tikhonov 1996; Sohn & Dvidge 1996; see also 
Hendry et al. 2005 for the latest review). Takada- 
Hidai, Aoki, & Zhao (2002) measured the total 
extinction toward the SN with a high-resolution 
spectrum and concluded that E(B — V) = 0.09 ± 
0.01. We adopt these values in this paper. 

Early-time optical spectroscopy and photome- 
try was published by Mazzali et al. (2002), Gal- 
Yam, Ofek, & Shemmer (2002), Kinugasa et al. 
(2002), Foley et al. (2003), Pandey et al. (2003a), 
Vinko et al. (2004), and others. Most authors em- 
phasized the spectral similarity to the broad-lined 
SN 1998bw and SN 1997ef. The early-time light 
curve (LC), was much broader than that of the 
normal Type Ic SN 19941, but was significantly 
fainter than that of SN 1998bw and peaked ear- 
lier. 

Yoshii et al. (2003) (hereafter Paper I) pre- 
sented early-time near-infrared (NIR) and optical 
photometry, obtained with the MAGNUM tele- 
scope. Using the UBVRIJHK photometry, they 
built a well-sampled OIR bolometric LC for SN 
2002ap. NIR photometry has also been reported 
by Nishihara et al. (2002) and Mattila, Mcikle, & 
Chambers (2002). Gerardy et al. (2004) described 
the evolution of NIR spectra, which were domi- 
nated by lines of intermediate-mass elements (see 
also Danziger et al. 2002). 

Mazzali et al. (2002) modelled the early-time 
optical spectra and bolometric LC using 1-D 
codes. They concluded that SN 2002ap was the 
energetic explosion of a C+O star, having evolved 
from a star of 20-25 M Q on the main sequence, 
and that it ejected ~ 2.5 M Q material with a ki- 
netic energy of ~ 4 x 10 51 ergs, including ~ 0.07 
M Q 56 Ni. This places SN 2002ap at the low-mass 
low-energy end of hypernovae. They also con- 
strained the explosion date to MJD 52300.0 ± 0.5 
days. We use the explosion date as the reference 
point for any SN epoch in this paper. Based on 
observed limits in pre-discovery images, Smartt et 
al. (2002) argued that the progenitor was either a 
single W-R star of main-sequence mass < 40M© 



or more likely part of an interacting binary. 

SN 2002ap was also monitored in the radio 
(Bcrgcr, Kulkarni, & Chevalier 2002; Soderberg, 
Nakar, & Kulkarni 2005), X-rays (Sutaria et al. 
2003; Soria, Pian, & Mazzali 2004), and UV (So- 
ria & Kong 2002). Berger, Kulkarni, & Chevalier 
(2002) estimated quite a small energy in any rel- 
ativistic material from the faintness of the radio. 
They addressed the material with velocities larger 
than the shock velocity (~ 0.3c). Therefore, their 
results do not contradict the high energy obtained 
in the 1-D model for the bulk SN ejecta by Mazzali 
et al. (2002), whose ejecta velocity distribution has 
a cut-off at 70,000 km s _1 . As a possible confir- 
mation of this result, the shock velocity derived by 
Bjornsson & Fransson (2004) from models of the 
radio and X-ray observations was also ~ 70,000 
km s -1 . 

SN 2002ap was not apparently associated with 
a GRB (Hurley et al. 2002). However, Kawabata 
et al. (2002) discovered a polarized spectrum com- 
ponent that resembled the total flux spectrum but 
redshifted by z ~ 0.3 (see also Leonard et al. 
2002; Wang et al. 2003). They proposed that this 
component is due to scattering from jet material. 
Totani (2003) further suggested that the jet may 
be too baryon-contaminated to trigger a GRB and 
that it was expanding almost freely and hence was 
radio quiet. Tominaga et al. (2004) reproduced a 
bump in the LC of the XRF 030723 afterglow with 
a SN 2002ap-like LC. 

SN 2002ap was monitored into the late nebular 
phase. Foley et al. (2003) published optical pho- 
tometry covering from 2002 June to December and 
spectroscopy continuing to 2003 February. The 
spectra show unusually strong and sharp-peaked 
[O I] and Mg I] emission lines (see also Leonard 
et al. 2002). Late-time spectroscopy was also 
obtained by other observers (Y. Qiu; K. Kawa- 
bata; S. B. Pandey; private communications), and 
late-time photometry by Vinko et al. (2004) and 
Doroshenko et al. (2003). 

It is also important to have intensive late-time 
NIR observations of Type Ib/c SNe. As demon- 
strated in the case of SN 1998bw (Patat et al. 
2001; Sollerman et al. 2002), the NIR contribu- 
tion to the total flux is quite large, in particular at 
late times. Late-time NIR observations may also 
reveal molecule and dust formation, if any exists 
(e.g., Gerardy et al. 2004). However, few Type 
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Ib/c SNe have been observed at late phases in the 
NIR. For SN 2002ap, the only late-time NIR ob- 
servation published so far is a spectrum taken in 
2002 August (Gerardy et al. 2004). 

In this paper, we report on BVRIJHK imag- 
ing photometry obtained for SN 2002ap between 
2002 June 12 and 2003 August 29, completing the 
1.5 years follow-up project with the MAGNUM 
telescope (see also Paper I) . We describe the obser- 
vations and data reduction and present our mul- 
ticolor LCs in §2. An OIR bolometric LC is then 
constructed using the photometry In §3, we dis- 
cuss the bolometric LC, analyze the evolution of 
the spectral energy distribution (SED), compare 
the late-time NIR LCs with those of other Type 
Ib/c SNe, and present 1-D models for the bolo- 
metric LC. 

2. Imaging Photometry 

2.1. Observations and Data Reduction 

Optical (Johnson-Cousins BVRI) and NIR 
(JHK) imaging photometry of SN 2002ap (012000 = 
01 h 36 m 23. s 85, (5 20 oo = +15°45'13."2) was car- 
ried out using the MAGNUM 2m telescope at 
the summit of Haleakala on the Island of Maui, 
Hawaii (Yoshii 2002). The multicolor imaging 
photometer mounted at the bent-Cassegrain focus 
has an SITc 1024 x 1024 pixel CCD with a scale 
of 0."277 pixel , for which the light is effectively 
received in a field of view of 430 x 430 pixel CCD 
or 119" x 119", and also has an SBRC 256 x 256 
pixel InSb with a scale of 0."346 pixel -1 yielding 
a field of view of 88."5 x 88. "5. A dichroic beam 
splitter enables simultaneous imaging through op- 
tical (UBVRI) and NIR (ZJHK'KL') filters 
(Kobayashi et al. 1998a). 

The observations were scheduled by the MAG- 
NUM scheduler (Kobayashi et al. 1998b), except 
for the solar conjunction from March to June in 
2003. Telescope dithering was performed with 6 
or 10-arcsec steps primary for JHK. The typical 
exposure time for one step was 190 sec for B, 95 
sec for V, 65 sec for each of RI, 60 sec for J and 
30 sec for each of HK. As SN 2002ap faded, the 
number of dithering steps was increased from 6 to 
12. The typical seeing size of the stellar image 
was 1.0-1.6 arcsec in V and 0.85-1.4 arcsec in K. 
Dome-flat images for BVRIJHK and NIR dark 
images were obtained at the end of each night. 



Image reduction was performed using our pack- 
age of IRAF-based 1 automated reduction soft- 
ware (MAGRED), which includes the standard 
corrections for bias and flat field frames. For more 
details see Paper I. 

Two reference stars, A and B, were observed 
in the same frame of SN 2002ap for differen- 
tial photometry. One reference star (A, 0:2000 = 
01 /l 36 m 19. s 49, (5 20 oo = +15°45'20."8) was used 
for BVRI as in Paper I. The other reference star 
(B, 0-2000 - 01 h 36 m 20. s 32, (5 2 ooo = +15°44'54."3) 
was used for JHK. In order to calibrate the 
magnitude of star B, NIR standard stars (Hunt 
et al. 1998) were also observed, whenever possi- 
ble, either before or after the observation of SN 
2002ap. The median magnitudes and the number 
of calibrations for star B in the InSb frame are 
J(6) = 15.69 ± 0.02, H(6) = 15.11 ± 0.05, and 
K(6) = 14.86 ±0.03. 

A variability check was made for stars A and 
B by increasing the total exposure time in pro- 
portion as SN 2002ap faded. Figure 1 shows the 
difference between their BVRI magnitudes after 
2003 July, when the optical flux of star B became 
detectable. (It was not possible to make the same 
plot in the NIR because star A was usually out- 
side or on the edge of the InSb detector.) The 
error bars estimated from photon statistics are in 
good agreement with the dispersion of the data 
points, so that it is not necessary to regard stars 
A and B as variable. 

Simple aperture photometry was used through- 
out, and the aperture magnitude of SN 2002ap 
was obtained as the magnitude difference between 
SN 2002ap and star A or B. An aperture size of 
6.9 arcsec was adopted for all images, which min- 
imizes the dispersion of the differential photome- 
try. However, special care was taken for reliable 
magnitude determination of SN 2002ap, given the 
presence of a faint background object within a 
distance of about 1 arcsec. This object became 
distinguishable from SN 2002ap, when SN 2002ap 
became sufficiently faint after solar conjunction in 
2003. In 2003 December, SN 2002ap faded com- 
pared to the background object, and it was easy 

X IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under con- 
tract to the National Science Foundation. 
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to resolve SN 2002ap and the object separately. 
Then, using the IRAF DAOPHOT package, the 
PSF magnitude of the object was measured sev- 
eral times when the seeing was relatively good. 
This way, the PSF magnitude of the object was 
determined as B = 21.7 ± 0.5, V = 21.6 ± 0.2, 
R = 21.25±0.12, / = 20.62±0.17, J = 20.0±0.1, 
H = 19.3 ± 0.1, and K = 19.1 ± 0.3. These fluxes 
were subtracted from each aperture photometry of 
SN 2002ap. 

Moreover, reference star B in NIR photome- 
try has a faint nearby star which is about 4 mag 
fainter. Photometry was performed on the images 
where this faint nearby star was subtracted using 
IRAF DAOPHOT package. 

2.2. Observed Light Curves 

The BVRIJHK data taken over 42 nights from 
June 12, 2002 to Aug 29, 2003 are tabulated in 
Table 1. The uncertainty in each magnitude was 
derived combining in quadrature the magnitude 
dispersion from the frames obtained during the 
same night and the magnitude error of reference 
star A for BVRI, or B for JHK. A small er- 
ror of 0.01 mag was added in quadrature for the 
BVRI magnitudes, taking into account the non- 
linearity correction error between SN 2002ap and 
reference star A. Another O.Olmag error caused by 
the seeing effect, arising from the profile difference 
between SN 2002ap and the reference star A, was 
also added. The absolute calibration error and the 
subtracted flux uncertainty of the background ob- 
ject were taken into account when evaluating the 
total photometry error. No color-term correction 
was applied to the magnitudes of SN 2002ap in the 
nebular phase, because strong line emissions and 
large color difference with respect to the main se- 
quence star make this correction particularly dif- 
ficult, as demonstrated by Menzies (1989). 

We show the MAGNUM U BVRIJHK LCs in 
Figure 2 {filled circles), combining the late-time 
data in Table 1 and the early-time data of Paper 
I. For comparison, we also show the UBVRI pho- 
tometry of Pandey et al. (2003a, b) (< 340 days; 
crosses), that of Foley et al. (2003) (< 320 days; 
open circles), and the early-time JHK photom- 
etry of Nishihara et al. (2002) (< 50 days; open 
circles). Obviously, our LCs have the largest time 
span. They are unique also in that both the early 
and late phase were well covered in NIR, the first 



time for any SN Ib/c. To guide the eye, we made 
B-spline fits (solid lines) to the data, excluding the 
Pandey and Foley late-time ones. 

It is clear from the figure that the decline of the 
late-time LC slowed down after September 2002 
in several bands. We derived the decline rates 
and the date when the decline rate changed in 
each band by x 2 -fitting the late-time MAGNUM 
data using a broken line (see Table 2). Before the 
change, which took place between about day 210 
and 250, the decline rate was similar in all bands, 
i.e., about 0.016 — 0.022 mag/day. The change is 
most profound in the NIR, as the J and H LCs 
flattened from 0.018 mag/day to 0.012 and 0.008 
mag/day, respectively. It is small for most optical 
bands, although the V LC did flatten from 0.021 
to 0.014 mag/day. As a result, the NIR contribu- 
tion to the total flux increased significantly going 
to very late times. This will be discussed further 
in the next section. 

Most of our late-time data agree well with the 
Foley et al. (2003) ones, although their B magni- 
tudes are 0.1-0.2 mag brighter than ours after day 
280 and their / magnitudes 0.1-0.25 mag brighter 
after day 160. The Pandey et al. (2003b) late- 
time data are 0.1-0.15 mag brighter than ours in 
the B band, 0.2-0.4 mag brighter in the V band, 
and 0.1-0.2 mag fainter in the / band. Those dif- 
ferences can be partly explained by measurement 
uncertainties inherent to practical broadband pho- 
tometry of SNe (e.g., see Hamuy et al. 1990), 
since these objects have time-evolving and line- 
dominated spectra of non-stellar nature. Our late- 
time BVRI LCs can differ from the those of Vinko 
et al. (2004) and Doroshenko et al. (2003) by as 
much as 0.3 mag. However, their data were less 
well sampled than ours. 

The late-time evolution of the optical colors, 
B — V, V — R, and R — I, is very different between 
SN 2002ap and SN 1998bw, as shown in Figure 
3. In particular, V — R of SN 2002ap between day 
140 and 330 is about 0.3-0.6 mag larger and B — V 
between day 200 and 340 about 0.3 mag smaller 
than in SN 1998bw. This can be easily understood 
within the context of spectroscopy. In the nebular 
phase, the spectra of both SNe were dominated 
by [O I] AA6300,6364 in the R band, [Fe II] mul- 
tiplets in the V and B bands, and Mg I] A4571 
in the B band, while Ca II] AA7291,7324 may also 
contribute to the R and / photometry and the 
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Ca II NIR triplet to the / photometry depending 
on the actual filter cut-off. As discussed by Foley 
et al. (2003), SN 2002ap has the strongest [O I] 
AA6300,6364 and Mg I] A4571 nebular emissions 
of any SN published. However, its [Fe II] multi- 
plets are much weaker than in SN 1998bw, a con- 
sequence of much less 56 Ni ejected by SN 2002ap. 
Between day 390 and 520, V—R and R—I dropped 
and B — V rose rapidly. This may indicate a flux 
shift from emissions of intermediate-mass elements 
to [Fe II] multiplcts. Unfortunately, no observed 
spectra have so far been obtained for Type Ib/c 
SNe at such late times. 

3. Discussion 

3.1. Bolometric Light Curve 

We built the OIR bolometric LC integrating 
the BVRIJHK broadband flux (see Paper I for 
details). A reddening of E(B - V) = 0.09 ± 0.01 
was corrected for (Takada-Hidai, Aoki, & Zhao 
2002), and a distance modulus of 29.5 was adopted 
(but sec Hendry et al. 2005). 

The LC is shown in Figure 4, extending that 
in Paper I to <~ 580 days since explosion. It con- 
sists of three independent data sets, which agree 
with one another. Most of the data {filled cir- 
cles), tabulated in Table 3 {late time) and Paper I 
{early time), were made from the MAGNUM pho- 
tometry. We interpolated and extrapolated the 
observations to estimate the the K magnitudes 
between day 280 and 340 and the J and I magni- 
tudes at day 520. For epochs < 11 days, the ob- 
servations of Mattila, Meikle, & Chambers (2002), 
Gal- Yam, Ofck, & Shemmcr (2002), and Cook et 
al. (2002) were used for interpolation. Bolometric 
magnitudes at day 407 and 578 {triangles) were 
converted from H magnitudes using approximate 
bolometric corrections. The other two data sets 
were based on the Foley et al. (2003) {open cir- 
cles) and Pandey et al. (2003a, b) {crosses) opti- 
cal photometry, respectively, which was combined 
with the Nishihara et al. (2002) and our JHK 
photometry. 

The bolometric LC at late time has a shape 
similar to that of SN 1998bw (Figure 4; squares), 
another hypernova, although the latter was much 
brighter. The SN 1998bw LC in Figure 4 is shifted 
down by 1.75 mag to match the peak magnitude 
of SN 2002ap. Both LCs show a similar late-time 



slowing-down. For SN 2002ap, the decline rate 
changed from ~ 0.018 mag/day between day 130 
and 230 to ~ 0.014 mag/day between day 270 and 
580. For SN 1998bw, the decline rate was ~ 0.019 
mag/day between day 70 and 220 and ~ 0.014 
mag/day between day 320 and 540. For compari- 
son, the late-time bolometric LC of another hyper- 
nova, SN 1997ef, followed the decay rate of 56 Co, 
i.e., ~ 0.01 mag/day (Mazzali et al. 2004). On 
the other hand, SN 2002ap reached its light max- 
imum ~ 5 days earlier than SN 1998bw did, al- 
though their peak widths are actually comparable 
(Figure 4; inset). 

3.2. Spectral Energy Distribution 

We show the SED of SN 2002ap at 5 typical 
late epochs in Figure 5, where the monochromatic 
fluxes converted from the MAGNUM BVRIJHK 
photometry are connected using spline-fitting 
curves. The observed upper limits of the if -band 
flux at day 393 and 520 are marked by arrows. 
Zero flux was assumed at both the the blue edge 
of the B band and the red edge of the K band. 
We used such SEDs to obtain the OIR bolometric 
magnitudes. 

The SEDs are dominated by the flux in the R 
and / bands before about day 340-390, which was 
mainly due to strong [O I] , Ca II] , and Ca II emis- 
sions (Foley et al. 2003). The B bump at the 
three intermediate epochs reflects the contribution 
of Mg I] A4571. After about day 340-390, the V 
flux, attributed to [Fe II] multiplets, increased rel- 
ative to that in other optical bands, while the I- 
band-dominated Ca II emissions died away. 

The NIR flux contribution rose rapidly after ~ 
day 300 (Figure 6). This can been seen in Figure 
5 as a big H bump in the SED and also in Figure 
2 as a significant flattening of the H and J LCs. 
To explain this, we examined the synthesized NIR 
spectra of the models that reproduce the late-time 
optical spectra of SN 2002ap (P. A. Mazzali et al. 
2005, in preparation). The most likely candidates 
are the strong [Si I] 1.646 fim and 1.608 (im lines, 
while the 1.099 /im line may account for the con- 
current J-band flux increase. The lines were al- 
ready developed in the day ~ 200 NIR spectrum 
of Gerardy et al. (2004), although they were still 
not as strong as a 1.5 ^m feature, possibly Mg I 
emission. We note in passing that the appearance 
of the CO first-overtone band in that spectrum 
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suggested by those authors seems to coincide with 
the -ftT-band LC steepening seen in Figure 2. 

The late-time iJ-band spectra of SN 1983N, 
taken by Graham et al. (1986) and of very low res- 
olution, showed a strong 1.65 /im feature, which 
was explained by those authors as [Fe II] lines but 
was soon, together with an adjacent feature, re- 
identified as [Si I] lines by Oliva (1987). These 
[Si I] NIR features can also be found in the syn- 
thetic late-time spectra of Type lb SNe computed 
by Fransson & Chevalier (1989). 

3.3. Late-time NIR Light Curves of Type 
Ib/c Super novae 

In Figure 7 we compare the late-time JHK LCs 
of SN 2002ap with those of other Type Ib/c SNe 
(see Table 4). SN 1998bw and SN 1984L were 1.5- 
2.5 mag brighter than all other SNe, indicating a 
large mass of 56 Ni ejected. The NIR LCs of SN 
1984L were relatively slow, as was its late-time 
optical LC which was modelled using very mas- 
sive ejecta but with normal kinetic energy (Swartz 
& Wheeler 1991; Baron, Young, & Branch 1993). 
The few points of SN 1998bw,when shifted down 
by 1.8 mag (J), 1.5 mag (H), and 1.7 mag (K), 
respectively, fall on the LCs of SN 2002ap. The 
LCs of SN 1983N have later and broader peaks 
than those of SN 2002ap, but they become simi- 
lar between in the J-band and if-band day ~ 200 
and 350. On the other hand, the JHK data of 
SN 19831 and 1982R seem to bridge the LC gap 
in SN 2002ap data between day 50 and day 140 
(solar conjunction). 

Figure 7 contradicts the picture that SNe 19831, 
1983N, and 1984L have similar JHK LCs (Elias 
ct al. 1985). Type Ib/c SNe were first established 
as a new class by those authors, who showed that 
these SNe have similar JHK LCs which arc dif- 
ferent from those of Type la SNe [and also by 
Wheeler & Levreault (1985) through an analysis 
of the optical spectra] . They artificially shifted the 
data of SN 19831 and SN 1984L in both phase and 
brightness in order to match the NIR LCs of SN 
1983N. Their SN 19831 data were assigned epochs 
20 days later than ours, and SN 1984L 14 days 
earlier. Our epoch estimates are more reasonable 
because we based them on optical LCs and spec- 
tra which are better observed and understood than 
the NIR data. 



3.4. Was SN 19831 similar to SN 2002ap? 

Inspecting by eye the only published spectrum 
of SN 19831 (Wheeler et al. 1987), we noticed 
that it resembles those of SN 2002ap between 
2002 February 16 and 22 (Kinugasa et al. 2002). 
The similarity can be seen not only in the overall 
shape and strongest line features, but also in rel- 
atively weak ones, like those near 4900A, 5900A, 
and 6200A. The spectrum was taken on 1983 May 
17, 5 days after discovery. We estimate the epoch 
of that spectrum as 22 ± 3 days since explosion, 
through this comparison, and hence derive an ex- 
plosion date of 1983 April 25. 

The LCs of SN 19831 seem to match those of 
SN 2002ap not only in the NIR (Figure 7) but 
also in the optical (Figure 8), if our estimates of 
the explosion date and the distance modulus arc 
correct. We adopt 31.0 as the distance modulus 
to the host galaxy of SN 19831, NGC 4051 in the 
Ursa Major cluster (Tsvetkov 1985). This value is 
also close to the average of the distance modulus 
of the Ursa Major cluster (31.4, Tully & Pierce 
2000), and that corresponding to its radial velocity 
corrected for the Local Group infall onto Virgo 
(30.5, LEDA 2 ; Paturel et al. 1997) 

Therefore, SN 19831 could be a precedent of 
the unusual SN 2002ap, which lies at the low- 
mass low-energy end of hypernovae (Mazzali et al. 
2002). However, its data are insufficient for us to 
make a solid conclusion. 

3.5. Modelling the Late-time Light Curve 

The late-time bolometric LC of SN 2002ap de- 
clines more slowly than that of the 1-dimensional 
(1-D) model of Mazzali et al. (2002) that best re- 
produced the early-time spectra and LC (see Fig- 
ure 9; dotted line). This is expected, since the LC 
follows that of SN 1998bw closely in terms of the 
decline rate, as shown in Figure 4. The best-fitting 
1-D model for the early-time spectra and LC of 
SN 1998bw also fails to explain its late-time spec- 
trum and LC (Nakamura et al. 2001; Mazzali et 
al. 2001). To mimic the outcome of 2-dimensional 
jet-induced explosion calculations (Maeda et al. 
2002), which applies to SN 1998bw and may also 
apply to SN 2002ap, Maeda et al. (2003) intro- 
duced a dense core. That core absorbs 7-rays ef- 

2 http://leda.univ-lyonl.fr/ 
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ficiently at late time. Using this structure, the 
slow LC decline of hypernovae was reproduced by 
those authors using a Monte-Carlo radiative trans- 
fer code but with simplified physics. 

We tested the dense-core scenario on SN 2002ap 
using our sophisticated 1-D SN radiation hydrody- 
namical and 7-ray transfer code (Iwamoto et al. 
2000; Mazzali et al. 2002; Deng et al. 2005). The 
best-fitting model (solid line) has ~ 0.6M Q ejecta 
below 3,000 km s _1 , compared with 0.1 M© in 
the model without a dense core, and ~ O.OIMq 
low-velocity 56 Ni. Compared with Mazzali et al. 
(2002), who modelled only the early-time obser- 
vations, the total ejecta mass has increased from 
2.5 M Q to 3 Mq, but with little change in the 
total 56 Ni mass (0.08 Mq) and the kinetic energy 
(4 x 10 51 ergs). A dense ejecta core is also required 
in 1-D late-time spectrum models (P. A. Mazzali 
et al. 2006, in preparation) in order to explain 
the observed sharp line cores of [O I] and Mg I] 
(Leonard et al. 2002; Foley et al. 2003). 

A dense core cannot be formed in 1-D explo- 
sion simulations for hypernovae (Nakamura et al. 
2001), but it is a natural product of 2-D jet- 
induced explosions (Maeda et al. 2002). This 
strongly indicates asymmetry in the SN 2002ap 
explosion, an intrinsic feature also shared by other 
hypernovae. Evidence of asymmetry in SN 2002ap 
can also be found in spcctropolarimetry (Kawa- 
bata et al. 2002; Leonard et al. 2002; Wang et al. 
2003), which shows an intrinsic continuum polar- 
ization varying around 0.5%, corresponding to an 
asphericity of ~ 10% for the bulk of the ejecta. 

We thank S. B. Pandey, E. Nishihara, and R. J. 
Foley for the electronic files of their photometric 
data. This work has been supported in part by the 
Grant-in-Aid for Scientific Research (16540229, 
17030005, & 17033002 for KN) and the 21st Cen- 
tury COE program (QUEST) from the JSPS and 
MEXT in Japan, and by the National Science 
Foundation under Grant No. PHY99-07949. 
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Fig. 1. — Difference between observed magnitudes 
of two reference stars [Am = m(A) — m(B)] in 
the BVRI bands for MAGNUM observations af- 
ter 2003 July. The mean difference (Am) is shifted 
to zero. Error bars are estimated from photon 
statistics. The standard deviation of data dis- 
tribution and mean error from photon statistics 
are 0.042/0.053 (B), 0.032/0.030 (V), 0.020/0.020 
(R), and 0.014/0.010 (/). The computed reduced 
chi-square values were 0.67 for B , 1.11 for V, 1.23 
for R and 1.33 for /. 

Fig. 2.— UBVRIJHK LCs of SN 2002ap for the 
first 1.5 years after discovery obtained with the 
MAGNUM telescope (filled circles). Also shown 
are the UBVRI data reported by Foley et al. 
(2003) (open circles) and Pandey et al. (2003a, b) 
(crosses) and the early-time JHK data by Nishi- 
hara et al. (2002) (also open circles). The LCs are 
shifted vertically in order to avoid overlap. Solid 
lines are B-spline fittings to the data, excluding 
the Pandey and Foley late-time ones. 

Fig. 3. — Evolution of optical colors, B — V, V—R, 
and R — /, of SN 2002ap compared with that of 
SN 1998bw. The colors of SN 2002ap are calcu- 
lated from the MAGNUM photometry, corrected 
for E(B -V) = 0.09, while those of SN 1998bw 
from the photometry in Galama et al. (1998), 
Patat et al. (2001), and Sollerman et al. (2002), 
corrected for A(V) = 0.1. 

Fig. 4.— OIR bolometric LC of SN 2002ap con- 
structed using the MAGNUM UBVRIJHK pho- 
tometry, assuming m — M = 29.5 and E(B — V) = 
0.09 ± 0.01, compared with that of SN 1998bw 
[z = 0.0085, m-M = 32.76, and A(V) = 0.1±0.1]. 
Also shown are bolometric magnitudes based on 
the Foley et al. (2003) and Pandey et al. (2003a,b) 
optical photometry, which are combined with the 
MAGNUM and Nishihara et al. (2002) JHK pho- 
tometry, and tho ones estimated from the MAG- 
NUM H photometry using approximate bolomet- 
ric corrections. The SN 1998bw bolometric LC 
is shifted down by 1.75 mag to match the peak 
brightness of SN 2002ap. It is constructed using 
the photometry in Galama et al. (1998), Patat et 
al. (2001), and Sollerman et al. (2002). 

Fig. 5.— SEDs of SN 2002ap at 5 typical late 
epochs, constructed using the MAGNUM photom- 
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etry. The BVRIJHK monochromatic flux is con- 
nected using spline-fitting curves. Capped arrows 
represent flux upper limits in the K band. 

Fig. 6. — Contribution of NIR flux to the total 
OIR bolometric flux in SN 2002ap compared with 
that in SN 1998bw. 

Fig. 7. — Absolute J {top), H (middle), and K 
(bottome) LCs of the Type Ib/c SNe whose late- 
time NIR photometry are available. See Table 4 
for data references and adopted SN parameters. 

Fig. 8.— Absolute B and V LCs of the Type Ic SN 
19831 (Tsvetkov 1985), assuming m - M « 31.0 
and an explosion date of 1983 April 25 (sec Table 
4), compared with those of SN 2002ap. 

Fig. 9.— OIR bolometric LC of SN 2002ap com- 
pared with the best-fitting LCs of the 1-D ejecta 
model with a dense core at v < 3, 000 km s _1 (to- 
tal M j w 3 Mo), and that without a dense core 
(total M C j w 2.5 Mq). Both models have total 
kinetic energy of 4 x 10 51 ergs and 0.08 M Q 56 Ni. 
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Data of light curves of SN 2002ap 
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Feb 23.22 


52693.22 


20.99±0.42 


















clear 


Feb 24.23 


52694.23 




20.49±0.11 
















clear 


Fob 25.23 


52695.23 






19.60±0.06 














clear 


Feb 26.23 


52696.23 








19.50±0.12 












clear 


Mar 11.21 


52709.21 
















17.89±0.13 




clear 


Mar 12.22 


52710.22 










19 


16±0 


22 






clear 


Jun 27.62 


52817.62 






21.26±0.22 














clear 


Jul 3.60 


52823.60 


22.80±1.36 


















clear 


Jul 5.60 


52825.60 




21.75±0.31 












18.80±0.21 




clear 


Aug 15.49 


52866.49 








< 21.5 












thin cloud 


Aug 28.50 


52879.50 










< 21.2 


19.83±0.34 




clear 


Aug 29.50 


52880.50 


















< 19.9 


clear 



Table 2 

Decline rates of light curves of SN 2002ap in the nebular phase 





B band 


V band 


R band 


/ band 


J band 


H band 


K band 


Bolomctric 


Epoch range (MJD) 
Decline rate (mag/day) 


52438.6 - 52533.5 
0.017±0.0004 


52437.6 - 52533.5 
0.021±0.0002 


52439.6 - 52817.6 
0.016±0.0001 


52441.6 - 52696.2 
0.017±0.0001 


52439.6 - 52508.5 
0.018±0.001 


52438.6 - 52524.6 
0.018±0.001 


52437.6 - 52577.5 
0.022±0.001 


52438.5 - 52533.5 
0.018±0.0003 


Epoch range (MJD) 
Decline rate (mag/day) 


52574.5 - 52823.6 
0.015±0.002 


52574.5 - 52825.6 
0.014±0.001 






52522.6 - 52710.2 
0.012±0.001 


52533.5 - 52879.5 
0.008±0.0005 




52574.5 - 52879.5 
0.014±0.001 


Date of change (MJD) 


52545±5 


52535±1 






52514±4 


52532±4 




52558±3 



Table 3 

bolometric magnitudes of sn 2002ap 



MJD 


a 


Day 


b 


M bol 




52438. 


.5 


136 


,5 


-13.53±0.03 


23% 


52458. 


5 


156, 


,5 


-13.17±0.03 


22% 


52478. 


.5 


176, 


,5 


-12.81±0.03 


22% 


52493. 


5 


191, 


,5 


-12.52±0.03 


22% 


52504. 


.0 


202, 





-12.32±0.03 


21% 


52508. 


5 


206, 


,5 


-12.28±0.03 


23% 


52521. 


.0 


219, 


,0 


-12.04±0.05 


21% 


52525. 


.0 


223, 


,0 


-11.96±0.03 


21% 


52533. 


5 


231, 


,5 


-11.83±0.03 


22% 


52574. 


.5 


272, 


,5 


-11.15±0.05 


23% 


52577. 


.5 


275, 


,5 


-11.14±0.05 


24% 


52601. 





299, 


,0 


-10.74±0.07 


25% 


52609. 


.5 


307, 


,5 


-10.62±0.07 


26% 


52617. 


.5 


315, 


,5 


-10.50±0.07 


26% 


52626. 


.5 


324, 


5 


-10.37±0.09 


29% 


52638. 


.5 


336 


,5 


-10.23±0.09 


29% 


52695. 


.0 


393, 


.0 


_q 47+0-17 

y -+' -0.18 


39% 


52709. 


.0 


407, 


,0 


-9.3± 0.2 




52822. 


.0 


520, 


,0 


-8.06t°i° 


54% 


52879. 


.5 


577, 


,5 


-6.8± 0.5 





a With an accuracy of 0.5 days due to data 
interpolation/extrapolation. 

b With respect to the explosion date (MJD 
52302.0) 
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Table 4 

Type Ib/c supernovae with late-time NIR photometry 





SN 1982R 


SN 19831 


SN 1983N 


SN 1984L 


SN 1998bw 


Spectroscopic type 


Ib/c 


Ic 


lb 


lb 


Ic, hypernova 


Distance modulus a 


31.2 


31.0 


28.25 


31.5 


32.76 


Extinction A(V) b 


0.5 


0.04 


0.2 


0.09 


0.1±0.1 


Explosion date (UT) 


1982 September 28 


1983 April 25 


1983 July 1 


1984 August 5 


1998 April 25.9 


JHK photometric 






Elias et al. (1985), 




Patat et al. (2001), 


data sources 


Mattila & Meikle (2001) 


Elias et al. (1985) 


Mattila & Meikle (2001) 


Elias et al. (1985) 


Sollerman et al. (2002) 



a /i = 0.72, fi A = 0.7, f! M =0.3; 1982R & 1984L, radial velocity corrected for LG infall onto Virgo (LEDA; Paturcl et al. 1997); 19831, see §3.4; 
1983N, Cephcids (Thim et al. 2003); 1998bw, z=0.0085 (Galama et al. 1998). 

b 1982R, Clocchiatti ct al. (1996); 19831, 1983N, & 1984L, Galactic Schlcgel, Finkbeiner, & Davis (1998); 1998bw, Patat et al. (2001). 

C 1982R, 26 days before discovery (Porter & Filippcnko 1987); 19831, see §3.4; 1983N, 2 days before discovery and 16 days before B maximum (Porter 
& Filippenko 1987); 1984L, 15 days before light maximum and by comparison with SN 1985F (Tsvetkov 1987; Schlcgel & Kirshner 1989); 1998bw, 
GRB 980425 (Galama et al. 1998). 



